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Abstract 

Using the scalar and tensor virial equations, the Lane-Emden equation expressing the hydrostatic 
equilibrium and small oscillations around the equilibrium, we show how the cosmological constant 
A affects various astrophysical quantities important for large matter conglomeration in the universe. 
Among others we examine the effect of A on the polytropic equation of state for spherically symmetric 
objects and find non-negligible results in certain realistic cases. We calculate the angular velocity 
for non-spherical oblate configurations which demonstrates a clear effect of A on high eccentricity 
objects. We show that for oblate as well as prolate ellipsoids the cosmological constant influences 
the critical mass and the temperature of the astrophysical object. These and other results show 
that the effect of A is large for flat astrophysical bodies. 

keywords: Large-scale structure of Universe. Galaxies: clusters: general. Instabilities. 

1 Introduction 

It is by now an established fact that the universe accelerates faster than previously anticipated El 
I^EIIS]. Hence some hitherto neglected ingredient (in general called Dark Energy) has to be responsible 
for this new phenomenon. To account for this phenomenon we can introduce new physics in terms of a 
scalar field [J] or modify three expressions in Einstein's equations which are often equivalent to a specific 
model with a scalar field. The three possibilities to account for the new physics are: the Einstein tensor, 
the energy-momentum tensor (this is to say, the energy momentum tensor of a fluid gets modified by 
the inclusion of other components, [H] or the equation of state [HI). The first possibility encompasses 
a positive cosmological constant A and higher order gravity with a more complicated Einstein-Hilbert 
action ^U]. In the the present work we choose to work with the cosmological constant as the simplest 
explanation for the acceleration of the universe. We shall put forward the question if such a cosmological 
model has an influence on astrophysical structures. We shall use equilibria concepts like hydrostatic 
equilibrium and virial equations to see how relatively low density astrophysical matter of different shapes 
behaves in a fast expanding universe. Anticipating the results, we can say that indeed there are some 
interesting effects. 

Often it is assumed that A does not have any effect on astrophysical processes which take place at scales 
different from the cosmological ones. Indeed, looking at the scales set by A, this assumption seems to be 
justified at the first glance. The scales set by the cosmological constant are of truly cosmological order 
of magnitude ^2^]). The density scale is set by A = SnGNPvac with pvac — (0.7 — 0.8)pcrit- The length 
scale, ta = I/a/A is of the order of the Hubble radius while the mass scale Ma = ta/Gn reaches up to 
the value of the mass of the universe. These scales constitute the so-called coincidence problem, namely 
the question as to why we should live exactly at an epoch where the scales of the cosmological constant 
are also the scales of the universe. Neither was it so in the past nor will it be so in the future when the 
universe expands further. The only astrophysical structures which match these scales are superclusters 
whose densities are indeed of the order of magnitude of pcrit- Indeed, here we can almost be sure that 
the cosmological constant is of relevance ^1 E). However, probing into astrophysical consequences 
of the cosmological constant of other, smaller and denser structures like clusters of galaxies or even 
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galaxies themselves, would look a hopeless undertaking unless we find circumstances where the effect 
of A (which in the very principle is present) gets enhanced. This can indeed happen through various 
mechanisms. For instance, in a problem where 7'a combines with a much smaller length scale, say ro, 
the effects can be sometimes expressed as rg r™. In consequence, the observable in which this expression 
enters gets affected by A in a way which is important at much smaller scales than r\. A concrete 
example is the Schwarzschild-de Sitter metric where we find the parameter ta together with the much 
smaller length scale of the Schwarzschild radius Tg. These two conspire in the form (r<,rA)^/'^ to define 
the largest extension of bound orbits as explained in the text and in |15| . We will discuss a very similar 
combination which emerges from the virial theorem defining the largest possible virialized structure 
with a given mass. Another possibility to enhance the effect of the cosmological constant is to consider 
non-spherical objects. It then often happens that the effect of A becomes (h/h)'' Pvac where k are two 
different length scales of a flattened object like a disk or an ellipsoid ^3]. This indeed happens for many 
astrophysical quantities, among other the critical mass, the angular velocity and the temperature (mean 
velocity of the components of the large structure) which we will discuss in the present paper. Finally, 
we can vary a dimensional variable to see if this enhances the effect of the cosmological constant. As 
an example let us quote the polytropic index n in the equation of state entering also the Lane-Emden 
equation. It is known that with growing n (n > 5) the object described by this equation of state does 
not have a well-defined radius as the density goes only asymptotically to zero. We will show that this 
pattern of behaviour becomes more dominant with A 7^ 0. 

Of course, all these effects become stronger the more diluted the mass conglomeration is. The superclus- 
ters are certainly the best candidates if we look for astrophysical effects of the cosmological constant. 
As a matter of fact, they do not seem to be virialized due to the extreme low density and their pancake 
structure jJB] where the effect of flatness mentioned above becomes powerful ^| E]. For the next 
structure, the clusters of galaxies (or groups of galaxies) with densities between one and three orders of 
magnitude above the critical density |17| . we would need one of the enhancing factors discussed above 
to see an appreciable effect of A. This is possible in various ways as shown below. Galaxy clusters can 
have various forms, among others oblate and prolate |18l I19| . And what is more, they can even rotate 
[TO] . We will show an explicit effect of A on their angular velocity and temperature in case the angular 
velocity is zero. Since the effect of large eccentricity is larger for prolate than for oblate ellipsoids, 
it is comforting to know that clusters can assume a prolate shape. For low-density galaxies like the 
Low Surface Brightness (LSB) galaxies whose density is roughly four orders of magnitude above the 
critical density |2UI 121) . we still find some effects. For ?? = 5 and A = the solution of the Lanc-Emden 
equation |221 is very often used as a phenomenologically valid description of the density profile (called 
also Plummer's law [23]. This is still possible as p ^ as r ^ 00. However, this property vanishes for 
low-density galaxies and the n = 5 case not only does not have a well-defined radius, its solution does 
not vanish asymptotically which thus rendering it unphysical. 

The paper is organized as follows. In the second section we will briefly review the general form of 
virial theorem including pressure, magnetic fields and, of course, the cosmological constant. Here we 
will also discuss some general results regarding A. In the third section we will specialize on spherical 
configurations. We will show how A sets the scale of a maximal virial radius and compare it to a 
result from the Schwarzschild-de Sitter metric. We will also solve the Lane-Emden equation numerically 
and analytically (for the polytropic index n = 2). In the fourth section we will discuss non-spherical 
configurations. First, we will show how A affects the angular velocity of spheroids. In addition, we will 
discuss the effects of A for the critical mass, mean velocity and mean rotational velocity too. The fifth 
section is devoted to small oscillations around equilibrium. 

2 Local dynamics with cosmological constant 
2.1 Newton-Hooke spacetime 

The cosmological constant enters the equations of Newtonian limit as a consequence of its appearance 
in the Einstein field equations. It is through this weak field limit approximation that A enters also in 
the equations describing the structure of astrophysical configurations. It is interesting to note that all 
variables to be found in the virial equations, are also present in the Poisson equation of the Newtonian 
limit. However, this is not always the reason why these terms enter the virial equations, at least in the 
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first order. The Poisson equation for a self-gravitating system modeled as an ideal fluid is written as 

V2<i> = 47rGAr(^p + 3^ + 2^^ -A. (1) 

where P is the pressure and Ucm is the electromagnetic energy density. The solution of at the zeroth 
order of v/c (from now on we set c = 1) is written as 

1>(r) = -G^ J^^ dV - lAIrp + • • • (2) 

where the dots stand for the correction terms that appear because the boundary conditions are now 
set at a finite distance jJJ. These terms can be usually neglected. The cosmological constant A > 
contributes to the expansion of the universe. This fact remains partly valid in the Newtonian limit 
where A gives us an external force. This defines the so-called (non-relativistic) Newton-Hooke spacetime 

[211221123 



2.2 The A-virial theorem 

The second order tensor virial equation can be derived in different ways: from a statistical point of view 
through the coUisionlcss Boltzmann equation, from a variational principle or by direct differentiation of 
the moment of inertia tensor 

^ik = / p{r)rirk d^r. (3) 



In the following we use the statistical approach |27] which also allows to derive higher order virial 
equations (for instance, the first order virial equation refers to the motion of the center of mass). In this 
context, from Boltzmann's equation one can derive the equation for momentum conservation (Euler's 
equation) written for a self gravitating system influenced by a magnetic field as 



+ p9.$ + ^5,(B2) + d, {V,j - B,B,) = 0, (4) 
where $ is the gravitational potential given by (0) (which includes A) and 

'Ptj = P{{vi - {vi)){vj - (vj))) = 6.,kP + TTik, Tr(7rjfc) = 0, (5) 

is the pressure tensor, P is the pressure and iTik its traccless part. Equations and Q together with 
an equation of state P = P{p, s) (s is the entropy) complete the description of a self gravitating fluid. 
By taking exterior products of r^ with Euler's equation and integrating over the volume of the system 
one obtains the second order virial equation as 

1 d^I, 



2 2T,fc-|Wfri + ^^GA.pvac2:.fc+n,fc where = J^P.kd\, (6) 

where Tik is the kinetic energy tensor and |Wf"^| is a generalized potential energy tensor which contains 
the contribution from the gravitational potential energy tensor Wj^. and the contributions of magnetic 
field through 

|WrNI>Vf,l(l-Aw). (7) 
The other quantities are defined as follows 

T^k^]: f p{v^){vk)d\, w!1 = -Gm f pir)rA^r)d\, A(,,) ^ (8) 

together with 

T,kiB) = 5,kB - 2B,k - f Tk l-6,,B^- - B,B, dS,, B,k I B,Bkd\, 6 = Tr(6,fc). (9) 

JdV J ^ Jv 
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A very useful version of the virial equation can be derived by assuming an isotropic pressure tensor and 
taking the trace in ©. This way we get the scalar A-virial equation 



1 d^X 



2 = 2^ - + gT^G^PvacX, (10) 



where the total kinetic energy is written as 

K.^- [ p(v^)d^r = T with U = [ Pd 

2Jv 2 ' Jy 



r 



(11) 



The equilibrium condition is reached for X = 0. This gives us the general A-virial theorem 

2T^k - iWfri + l^P^^cX^k + =0, 2/C - iW^cnj ^ ^^p^^^j = 0. (12) 

For rotating configurations with constant angular velocity, the kinetic term is modified as in the standard 
way as 

Tik ^ Tik + TZik, T^ifc = 2 (^rot^ifc ~ ^rotJfcj^^™') , R = Tr(7?,ij), (13) 

with TZij the rotational kinetic energy tensor and Tik is referred to motions observed from the rotating 
reference frame. The A-virial theorem has been used in different contexts in |12l 1131 f?51 129| . In the 
present work we will extend these studies. 



2.3 General consequences 

The tensor virial equation is widely used in many astrophysical applications. The inclusion of A provides 
a new way to study effects of the cosmological constant (parameters , in general) on astrophysical 
objects. The outcome depends essentially on two factors: the geometry of the configuration and the 
density profile. Wc will explore the spherical geometry for both constant and varying density profiles 
and study some effects for non spherical geometry with constant density. The consequences that can 
be derived from the A-virial theorem can be classified in two categories. The first one puts an upper 
bound on the cosmological constant or alternatively a lower bound on density of objects in gravitational 
equilibrium. Provided these bounds are satisfied, we can also study in the second step the effects of A 
on other properties of the astrophysical configurations like rotation, small oscillations etc. 
The first simple consequence of the virial equation emerges if we require the system to satisfy (|12ll . The 
fact that /C > implies an upper bound on the vacuum energy density 

. ^ 3 ^^^^ 



"•"^ - SttI Gn 

All systems in equilibrium have to satisfy (|14|) . Note that the right hand side of this expression is a 
function of both the density and the geometry of the system. Hence, we must expect different bounds 
for different geometries and density profiles. For instance, if we assume a constant density and B = 0, 
we can define VV^ and X through 

iw^l = ^gVIw^'I, T^pi, (15) 

such that the bound written in (|14() becomes 

P>-4pvac, with A=^i-^\. (16) 




The factor A which is only a function of the geometry (if we neglect the contribution of magnetic 
fields), will appear in many places in the paper. Its relevance lies in the fact that it enhances the effect 
of the cosmological constant for geometries far from spherical symmetry when A is large. A useful 
generalization can be done for situations in which we use the tensor form of the virial equation, namely 

16n I Xij \ 

" = — Uw^l ] ' ^ " constant. (17) 
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Finally, a curious equation can be derived by eliminating A from the tensor virial equations: 

2?;.-|nf"l + n., 

2T -iws™i-i-n J ' ^ ' 

Although A does not enter this equation, p8|l is only valid if the denominator is non-zero as is the case 
with A ^ 0. In Sect. 4 we will use this equality to infer a relation between the geometry and rotational 
velocity of an ellipsoid. Having discussed the general form of the virial theorem, we will discuss now the 
effects of A and set B = 0. 



3 Spherical configurations 



The tensor virial equation is trivially satisfied for spherically symmetric configuration without a magnetic 

ik 



field, since — (5ifcVV^ and lik = SikT. Therefore, in this section wc use only the scalar form of JB)). 



3.1 Constant density 

Explicit expressions can be derived in the spherical case with constant density, with jW 



Nl _ 3 GnM 



' 5 R 

and X = |Mi?^, so that -4sphcricai = 2. In this case the ratio p^acj P does not get enhanced much by the 
geometrical factor A. 

Worth mentioning is the result from general relativity. There the upper bound for the cosmological 
constant comes out as A < AttGnP |3UI where p is the mean density defined hy p = M /V . This 
bound is derived not only from Newtonian astrophysics, but also from a general relativistic context via 
the Tolmann-Oppenheimer-Volkoff equation [32] for hydrostatic equilibrium of compact objects |12II31| . 
Another relevant effect of the cosmological constant is the existence of a maximal virial radius of a 
spherical configuration which can be calculated from the A-virial equation. Using the expressions for 
|yy^| and X given before, equation H10|l yields as a cubic equation for the virial radius i?vir 

Rl„ + {iQwD fivir - ir.rl = 0. (19) 

Here we introduced the dimensionless temperature rj as 

— = T, (20) 

where /x is the mass of the average member of the configuration, /cg is the Boltzmann constant, T is the 
temperature and is defined by 

r, = GnM, (21) 

The length scale ta is set by the the cosmological constant as 



TA = = 2.4 X 10^ r!-i/2 Mpc « 1 



X 10^° ly. (22) 



hjQ « 0.7 is dimensionless Hubble parameter and l^vac = Pvnc/ Pa-it ~ 0.7 is the density parameter 
at the present time. The positive real root of equation p9(l is given by 

-Rvh(r/) =Tn(77)i?vir(0), (23) 

where i?vir(0), is radius for the configuration at 77 = is given by 

i?™.(0) = {3xf^ TA = (3r,rA)l/^ (24) 
and the dimensionless parameter x is defined as 

a; ^ II = 1.94 X 10-23 (1L\ hj^nlil « 1. (25) 
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The radius i?vir (0) is the largest radius that a spherical homogeneous cloud may have in virial equilibrium 
(i.e, satisfying H12|l '). The function ■07(77) can be obtained from the solution of the cubic equation and 
reads (clarified in the appendix A) 



Tu{7j) = 2.53a;-^/^7//2 sinh 



iarcsinh (0.24x1] 



(26) 



Figure 1 shows the behavior of tn(r;) for different values of x. We sec that the increase of the temperature 
implies a decrease of the effects of A which can be easily checked if we solve i?vir from the virial theorem 
with A = and compare it to the approximation 77 ^ 00 in H23|l : 

ii-vh.(A ^ 0) = R^Uv ^ c») = R^M = — • (27) 

10 77* 

We can consider (|23() as a radius-temperature relation for a fixed mass applied on astrophysical structures 
in a single state of equilibrium in the presence of A. That is, given x and rj we calculate the radius. But 
we can adopt another point of view for this relation. Imagine a spherical configuration characterized by 
a constant mass M . In analogy to a thermodynamical reversible process, the configuration may pass 
from one state of virial equilibrium to another following the curve vo — rj, that is, satisfying the condition 
X ~ Q. Clearly, there must be some final temperature 7/^, when this process ends since the temperature 
cannot increase indefinitely. But of course since the virial equations are not dynamical we cannot know 
which stage is the final one. If we assume that the effects of A are negligible when rj ^ rj^,, then using 
Eqs I^nil and ^ we get 

77. ^ 0.208^7;lx2/^ = (28) 

This is an equation for the temperature T^, 

M-2/3j,^ = 0.138 nj-^rl^^^. (29) 

For a hydrogen cloud (/i ~ yTiproton); we then write the mass-temperature relation using Eas. (|22f) and 
(1211) as 

/ n,r \ 2/3 

= 8.60 X 10-4^-1 — h%'nli^ K. (30) 



Note that this expression maintains the same dependence of the standard mass-temperature relation 
derived from the virial theorem, i.e, T oc M^/'^ (see [221 or equation H71|l of this paper). However the 
meaning of H30|l is different from that of typical mass temperature relations since (|3()() is associated 
to the temperature that a system acquires in the final stage after after going through some reversible 
processes which took the system through successive states of virial equilibrium with constant mass from 
a radius i?vir(0) to a radius i?vir('7*) or vice versa. On the other hand, the mass-temperature relation 
like Eq. (|71() of this paper relates the temperature of any configuration in equilibrium with its observed 
mass at constant density. In that context one considers only one equilibrium state and the cosmological 
constant enters just as corrections. 

As a final remark on Eq. (|23(l . we discuss a result which formally coincides with the virial radius i?vir(0) 
derived from the Schwarszchild de-Sitter spacetime but whose physical meaning is quite different. 
The Schwarzschild-de Sitter metric takes the form 

= - (1 - ^ - ^) di' + (1 - ^ - ^) ' d7'' + '-'d^^' + sin' ^d^'> (31) 
Now in contrast to H21|l we have 

Rs = GN/i, (32) 

with /.J the mass of the object giving rise to the Schwarzschild-de Sitter metric (in Eq. (|21(l M is the mass 
of the total conglomeration whereas here we consider fj, as the mass of its average member). Choosing 
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Figure 1: Ratio between -f?vir('?) and Rvir(O) for different values of a; as a function of r; = {3/fi)T, where fj, is the mass 
of the main average components of the system. 



the afRne parameter as the proper thne t the equation of motion of a test-body can be cast in a form 
similar to the corresponding equations from non-relativistic classical mechanics. 

1 / dr \ ^ 1 

2 ( d7 ) +Uctf^^{£^-l)^C^ constant, (33) 

where £ is a conserved quantity and [/off is the effective potential, defined by 

, ss di , , Rs 1 , , 

£ = {1 + 2[/eff(r)) — , C/eff(r) = ^ - « — • (34) 

dr r 

For simplicity we are have chosen here the angular momentum L to be zero. With L zero or not, Ucs 
displays a local maximum below zero due to A 7^ forming a potential barrier. This is to say, the 
standard local minimum where we find all the bounded orbits is now followed by a local maximum after 
which Uoff goes to —00. With A = this function approaches zero asymptotically. One is immediately 
tempted to say that this barrier will occur at cosmological distances. This is not the case and one 
calculates 

W = {SRsrlf = 9.5 X 10-^ (t7-X^' ( P^Y^\-y-' Mpc. (35) 

V^WQ/ VP vac/ 

In other words, the combination of the large scale r\ with the small scale Rg gives us a distance of 
astrophysical relevance, namely Tmax- Its relevance lies in te fact that beyond Tmax there are no bound 
orbits. Indeed, with fj. the solar mass, rmax is of the order of a globular cluster extension (70 pc); with 
II as the mass of globular cluster, r^ax comes out to be of the order of the size of a galaxy (10 kpc), 
and finally taking /x to be the mass of a galaxy, rmax gives the right length scale of a galaxy cluster (1 
Mpc). Certainly, the value of the extension of a large astrophysical body is the result of a multi-body 
interaction. But with the actual values of Tmax, it appears as if the length scale (we emphasize that 
we are concerned hare about scales and not precise numbers) of an astrophysical conglomeration is 
approximately rmax, which apparently means that this scale does not change drastrically when going 
from a two body problem to a multi-body calculation. This makes sense if the object under consideration 
is not too dense. We can now say that whereas M in i?vir(0) (via rg) is the mass of the object, /i in Eq. 
(|33|) is the mass of its members. Clearly, we have i?vir(0) ^ J'max, but both scales are of astrophysical 
order of magnitude. A result related to (|24|) derived in the framework of general relativity can be found 
in inj. 

3.2 Non-constant density 

The examination of configurations with non-constant densities can be done in two directions. Knowing 
the density profile p{r), we can set up the virial equation and evaluate the equilibrium conditions from 
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the inequality 114|l . In this picture, the effects of A are included in the solution for the potential $ as 
in Eq. |(2J) and the resulting term acts like an external force, as mentioned before. 

The second option is to combine the Eqs. 1^, Q and an equation of state (e.o.s) for which we can take 
a polytropic form P = Kp^^^ . This way we obtain the Lane-Emden equation with A ^2] 



1 d 



>d^ 



where pc is the central density, r 
Jeans length defined as 



Pvac 
Pc 

a^, p{r) = peV'"(C) with V(0) 



lK{n + 1) 



Anpc 



(36) 

1, ip'{0) =0 and a is the associated 



(37) 



It is important to notice that in this picture the expected effects of A are to be found in the behavior 
of the density profile since now Eq. (|36() implies that its solution is also function of the parameter C.c- 
Some effects of A are contained in the total mass and the radius of the configuration which is reached 
when ip{^i) ~ 0. Hence (|36ll implies 



df V dC 



-I 1/2 



R = a^i = VA 



K{n + l)pc" 



de d^J^, 



1 1/2 



(38) 



Note that the radius of the configuration is now proportional to r\. This is due to the fact that A 
sets a scale for length. However, this does not mean that R will be always of the order ta as A is also 
contained in the expression in the square brackets in Eq. H38|l . We will show this below in a concrete 
example. Since A is a new constant scale the Lane-Emden equation loses some of its scaling properties 
as explained in |12) . The mass of the configuration can be determined as usual with, 



Af(C) 



Ana^Pc 



Cc 



3di/; 
?d?. 



(39) 



where we used Eq. H36|l for the second equality. The total mass is then obtained by evaluating the last 
expression at ^ = ^i. As expected, the mass increases because the Newtonian gravity has to be stronger 
in order for the configuration to be in equilibrium with A 7^ 0. Figure 2 shows the numerical solutions 
for n = 1 to 71 = 5. We expect that the radius of the configuration is increased by the contribution 
of Cc and find it confirmed in the figures. However, not always is the radius of the configuration well 
defined, even if n < 5. For sizeable values of Cc (black line) we cannot find physical solutions of Eq. 
()36() as the function ip acquires a positive slope. One might be tempted to claim that the radius of the 
configuration could be defined in these situations as the position where ip has its first minimum, but as 
can be seen for n = 3 such a radius would be smaller than the radius with Cc ^ which contradicts the 
behaviour shown for the other solutions where ^(Cc 7^ 0) > ^(Cc = 0). As already mentioned above this 
is the correct hierarchy between the radii because large Cc gives rise to a large external force pulling at 
the matter. The numerical solutions show that for relatively large values of Cc only n = I has a well 
defined radius. In this case the effect of A is a 13% increase of the matter extension as compared to 
Cc = 0. As we increase the polytropic index, Cc ~ 10^^ leads to non-physical solutions while the effect 
with bigger values of Cc becomes visible only for ri = 3. For instance, Cc ~ 10~^ results in a radius 
which is 17% bigger than the corresponding value with Cc ^ 0. The combination n ~ A and Cc ~ 10"'^ 
also leads to a non-physical solution, whereas the radius of the case Cc ~ 10"'' displays a difference of 
13% as compared to Cc 0. Finally, for n = 5, the only physical solutions are obtained for the lowest 
values of Cc where ip' < This case is particularly interesting as with A = it is often used a a viable 
phenomenological parametrization of densities |23l I35| . The solution has an asymptotic behaviour as 
r~° which has been also found in LSB galaxies de Blok at al. 2004). With A 7^ the n — h seems less 



appealing as the matter is diluted. For all values of n, 
negligible. 



the difference between Cc = 10 and 10 is 



Analytical solutions of (|36|) can be found for 71 = 0, 1 and 7i = 5 if A = 0. For n 00, the polytropic 
equation of state reduces to the equation of state of the isothermal sphere P = up. As an example, for 
A 7^ 0, we can write the analytical solution in the case n = 1 as 



Cc)^+Cc. 



(40) 
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Figure 2: Effects of A on the behaviour of the density of a polytropic configuration for different ratios and different 
polytropic indices. The radius of the configuration is not always defined, even for n < 5. For higher values of pvac, only 
the n = 1 case has a definite radius for these values of fc- For other cases, the configuration is defined only for small C,c- 



The radius is i? = a^i, where is the solution of the transcendental equation 

CiCc = -(1 - Cc)sin^i. 
In the first order of Cc one finds 



R 



Tin (1 + Cc) . 



(41) 



(42) 



Equation l|41(l also implies that there exists some Ccrit such that for C,^ > Ccrit, we cannot find a real 
solution for . Approximately this gives 



Pc > 10. 



(43) 



which, provided the overall density is not too big, is better that p > 2pvac which is a result from the 
general inequality H16(l for p = const and spherical symmetry. Finally, we can calculate the contribution 
of A to the total energy of the object. Generalizing the results found in we obtain 



6 it V 5 — 71 



R 



(44) 



from which we infer that the correction is very small in this case. 

At the end of this section we would like to summarize the findings from Fig. 2. In table 1 we write the 
ratio = 0)/^i(A ^ 0) for the same ratios Cc and polytropic index as in Fig. 2. The horizontal line 
represents a non-defined radius. The symbol oo indicates that the radius is defined only asymptotically 
in case of A = 
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Cc 


n = 1 


n = 3/2 


n = 2 


n:=:3 n = 4 n = 5 


0.1 


~ 0.88 








0.005 


~ 1 


~ 1 


- 0.98 


- 0.86 


2 X 10-4 


- 1 


- 1 


- 1 


- 1 - 0.88 


1.2 X 10^5 


- 1 


- 1 


- 1 


^1 ^1 


6.1 X 10-^ 


- 1 


- 1 


- 1 


- 1 - 1 oo 



Table 1: Values of the fraction 5i(A = 0)/^i(A ^ 0) for different values of the ratio and the polytropic index. The 
horizontal lines represents the non well defined radius. 



4 Nonspherical configurations 

4.1 Rotating configurations 

As emphasized already before, the effect of A can get enhanced for non-spherical objects. This happens 
when the vacuum energy gets multiplied by a ratio of two length scales li and I2 and we end up 
with expressions like Pvacih/h)" ■ For instance, for ellipsoidal configurations the geometrical parameter 
A entering among others the inequality H16|l can be calculated from its definition H17|l with = 
-^T:aia2a3Sija^ and Wik given in We have jj^l 



A, 



obi 



4 e(3 



arcsm e 
2e2 



3(l_g2)3/2 



2yr 



In 



37r 



1 



Oi 

as 

aa^ai 2 

^ 3 



In 



203 
ai 



(45) 



where the eccentricity^ e is = 1 — a|/a^ for the oblate and = 1 — af/a'^ for the prolate case. The 
behaviour of .Aobi and ^pro is shown in figure UTI It is clear that for very flat astrophysical objects 
we can gain in this way several orders of magnitude of enhancement of the effect of A if ^ is a factor 
of attached to pvac. Needless to say that we often encounter in the universe flat objects like elliptical 
galaxies, spiral disk galaxies, clusters of galaxies of different forms and finally superclusters which can 
have the forms of pancakes. Of course, the more dilute the system is, the bigger the effect of A. We can 
expect sizeable effects for clusters and superclusters, even for very flat galaxies. Regarding the latter, 
low density galaxies like the nearly invisible galaxies are among other the best candidates. 
A convenient way to model almost all flat shaped objects is to consider ellipsoids which in the limit of 
flattened spheroids can be considered as disks. There are three different kinds of elliptical configurations, 
characterized by three semi-axes ai = a, a2 ~ b and a^, = c: oblate , with a = b < c, prolate with 
a = b > c and triaxial systems with a > b > c. Here the tensor virial equation provides a tool to 
determine which of these geometries are compatible with the virial equilibrium. Considering the case 
A = and p =constant, or spheroids with confocal density distribution whose isodensity surfaces are 
similar concentric ellipsoids |35II36| i.e. 



p — p{m'^), rri^ 



3 2 

xf 



(46) 



the oblate ellipsoid (MacLaurin) emerges as a solution of the virial equations with a bifurcation point to 
a triaxial Jacobi ellipsoid |37] . We can view the virial equation without A as a homogeneous equation. 
Switching on A ^ this becomes an inhomogcncous equation whose right hand side is proportional to 
Pvac- It is therefore a priori not clear if in the case of A 7^ we can draw the same conclusions as with 
A = 0. 

Let us assume a configuration which is rotating with constant angular velocity around the z axis. 
Neglecting the internal motions, the A-tensor virial Eq. for such a configuration is 



^rot i^ik ~ 5izlz 



ik\ 



-TrGNPvac^ik = -(^ifcll. 



(47) 



^Once the density is given, the inequality 1161 becomes in this way a defining equality for Cmax such that e < Cmax in 
order to maintain equilibrium. 
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Figure 3: Function g{e) and geometrical factor ^(e) for prolate and oblate ellipsoids. These functions have their largest 
values for flat oblates and large prelates. 

The equations with diagonal elements yield the following identities 



1^1 



-ttGnPv 



$1^ T 



This set of equations can be resolved for the angular velocity: 



- 

"rot — 



-TrGArpvac ( 1 - 



The same expression holds if we make the replacement 



Wyy and Jx 



Wzz\+o^Pv.cIzz- (48) 



(49) 

Xyy on the right hand 



side of H49|l . We can also eliminate ri^^^. from these expressions to arrive at the condition 



|w.?J-|w. 



-GnPv 



(50) 



with 5 1 = Tyy/Xxx- This expression determines the possible ellipsoidal configurations. One solution is 
obvious: (5/ = 1 together with IW^^.! = This happens if the density is constant or for a confocal 

ellipsoid with the density given in H46|l . Hence this is not different from the A = case. However, the 
bifurcation point to a triaxial system will get affected by A Furthermore, the explicit calculation 
of this bifurcation point will now depend also on the details of the density profile even if we take the 
latter to be as in Eq. This is a direct consequence of a 'pre-existing' density scale pvac- For a 

homogeneous configuration we have 



- _ /l-i (l-ri ) 



1 - 



Pvac 

p 



(51) 




where Sj = lyy/l^x, Vi = |W^|/|yV,j^ | and An arc defined in Eq. lfT7|l . For an oblate configuration they 
can be calculated in terms of the eccentricity e = ^1 — a? j (? to give |35| 



Axx — 



VV 



3^/1 



arcsm e 



3^/^ 



arcsm e 



Furthermore we have 

r]x^'ny= 2(1 - 6^)1/2 



arcsm e 



(52) 



(53) 
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Using these expressions, we write the constant angular velocity for an oblate ellipsoid from the first line 
of equation Eq. as 



"rot — "O 



— 3(e) 



where Vlo corresponding to the angular velocity when A = is given by the Maclaurin formula 



-.■^■^xx (1 - Vx) = 27rGjvP 



^ „ ' (3 -2e2) arcsine - 4(1- e^) 



(54) 



(55) 



The function g{e) defined through Eq. (|51|) can also be calculated with explicit dependence on the 
eccentricity as 



.9(e) 



VxAz 



'II = Zg5 

I - Vx 3 



(1 - 6^)1/2(3 - 2e^) arcsine - 3e(l - e^ 



(56) 



As is evident from the above equations, A has a twofold effect on the angular velocity. Firstly, it reduces 
the angular velocity with respect to the value f2o especially at the local maximum (see Figure 3). This 
is not a small effect and can affect even galaxies. Secondly, we see from Eq. H55|l that fip ^ for e — > 1. 
On the other hand 



-.9(e) 



Pvac 

327r / 03 



9 V ai 



(57) 



approaching 1 for a very flat oblate configuration and not too dense matter. Therefore, beyond the 
local maximum in fij-ot the cosmological constant causes a steeper fall of firot towards 0. Another 
relevant interesting quantity which can be calculated in this context is the ratio of the rotational over 
the gravitational energy contributions to the scalar virial equations, i.e., 7?./|yV^|. In accordance with 
Eq. 1)55(1 the latter can be written as 



/3 = 7^/|W 



f3 = Po 



Pvac 

p 



.9(e) 



3 

2^ 



1 - 



e(l 



=211/2 



arcsm e 



- 1. 



(58) 



The effects on /3 are therefore similar to the the ones encountered in fJiot (see also Figures 3). Finally, 
on account of P > we can infer from the virial equations with p =const the following inequality 



which together with H58|l results in an inequality for the density of an ellipsoidal configuration: 



p > Apvac, 



^ 1 fA-4f3og\ 



(59) 



(60) 



For e ^ 0, we have A 2 A, while for e ^ 1, .4 ^ ^, and A> A. Therefore the above inequality is 
slightly weaker than the bound given in 1(16(1 . Nevertheless it is useful since it it derived directly from a 
different starting point (P > 0) than 1(16(1 which is based on /C > 0. 



4.1.1 Mean mass- weighted rotational velocity 

Not always the deviation from spherical symmetry guarantees that the effect of A on observables is 
sizeable. This depends on the context and also which scales we compare. If we compare pvac to densities, 
the quantities A and An defined in Eq. ((16(1 are for flattened objects large enough to enhance the effect 
of A. If r\ is combined with the Schwarzschild radius Ts to give rmax as in Eq. ((24() . the result is still of 
astrophysical relevance. But if we had to compare one of the axes of an ellipsoid to ta (i.e. ai/r\) the 
effect would be negligible unless the extension is bigger than Mpc (clusters and superclusters) and 
the small ratio a^/rA is comparable to another small quantity of the same order of magnitude entering 
the equation under consideration. This happens for instance if we generalize a result (discussed jSHI) on 
a mass-weighted mean-square rotation speed of an ellipsoidal object to include A. Assume that due 



12 



to symmetry properties of the object the only relevant components of the tensors in the tensor virial 
equations are xx and zz. We then obtain 



If the only motion is a rotation about the z axis we have Tzz = and we can solve for T^x as 



(61) 



1 



1 



(62) 



Using Eqs. (O. © and (jSJ, the quantities entering our equation can be parametrized in the following 
way 



2r,. = \j p{vl)d\ ^ ]^Mvl n,, - Mai ^zz = (1 - ^o)n. 



(63) 



where Vq is the mass weighted mean angular velocity, cto is the mass-weighted mean-square random 
velocity in the x direction and So measures the anisotropy in 11^^. If So is of the order of one, it suffices 
to compare ctq with AJ^x/'^M. If both are of the same order of magnitude, the effect of A is non- 
negligible. Since vq and cto are of the same order of magnitude and vq is non-relativistic, we can assume 
that (To < 10^^. The quantity AIxx/'iM can be estimated to be (ai/rA)^. Hence if ai ~ IMpc, ■^o (co) 
has to be truly non-relativistic and of the the order of 10~^ to gain an appreciable effect of A. This 
improves if is one order of magnitude bigger which is possible for large galaxy clusters. The velocities 
have to be then at most 10~'*. In these cases we have to keep A and while solving (I63f) for Uq one has 



1^ 
2a„^ 



il-5o) 



\W^x\ 



1 



1 Air 



3 Mcr^ 



I.. \W 



N I 

XX I 



^xx |WN I 



1 



(64) 



Note that if A cannot be neglected Vq / CTq is not only a function of the eccentricity for ellipsoids with the 
density give in (|46|) , but depends also on the details of the matter density as the latter does not cancel. 
Using almost the same set-up as above, we can use equation (|18|l to establish a relation betweenuo, the 
mass M and the geometry of the object which we choose below to be oblate. After straightforward 
algebra we obtain 



2ai 



3 Gat il/VT^ 



10 



3 arcsin e 



+ (1 - So) 



- 1. 



(65) 



Since this relation is derived from H18|l which in turn is based on the assumption of A 7^ it is only valid 
for non-zero cosmological constant albeit the latter does not enter the expression. Note the enhancement 
factor (ai/as)^. 



4.2 Other effects of A for non-rotating configurations 

We can now derive other relevant quantities from the scalar A-virial theorem applied to homogeneous 
ellipsoidal configurations. In this section we will not consider rotating configurations, but systems with 
kinetic energy coming from internal motions. We will focus again on ellipsoidal, oblate and prolate, 
geometries. As in the preceding section the relevant quantity here is the function A written for both 
configurations in equation (|45|) . 

Critical mass: Consider the stability criteria for a homogeneous cloud with mass M and internal mean 
velocity (v) in mechanical equilibrium with the background with pressure P (see j39p. The system will 
collapse under it's own gravity (P < 0) if its mass is greater than a critical mass Mc. With cosmological 
constant, this critical mass is increased with respect to its value Mcq when A = 0, which is expected, 
since now there is an external force associated to A that acts against Newtonian gravity and hence the 
collapse can be postponed. By using the scalar virial theorem (I12() and setting P = as the criterion 
for the onset of instability, we can write for arbitrary geometry 



Mr_ = M, 



cO 



I- A 



PcO 



1 -1 
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PcO = 



Mco 
V 



(66) 




Figure 4: Effects of A on the angular velocity f^j^t and the ratio /3 for pvac ~ O.Tpcrit 



This expression is valid for any geometry. However, we pointed out already that spherical symmetry 
implies A ~ 2 and the effect is suppressed. Some numerical values can be considered by writing 
(w^) = SksT/mp for a hydrogen cloud with T « 500 K and radius R « lOpc. One then has 

A ( ^) « 10-8, (67) 

V PcO J 

which represents a very tiny correction to the critical mass Mco for a spherically symmetric object. On 
the other hand, for ellipsoidal configurations with the same temperature, we have 

SM,A..,{S^)^...n"{^)\ .M = A.,(^)»10-.(j^)'. m 

in the e ^ 1 approximation for oblate and prolate configurations, respectively. We have set c = lOa for 
the prolate case and f^vac = 0.7. For an oblate ellipsoid with a « 50 kpc, one has SM « 0.15 while for 
the prolate ellipsoid with c « 50, SM « 10^. 

Mean velocity and Mass- Temperature relation. By virtue of the scalar A-virial theorem, we can also write 
down the mass-temperature relation for an astrophysical structure. Note first that using JC = ■^AI{v'^) 
in equation (|12|) we have for A = the standard expression for the mean velocity 

>^=° = ^ = M ' 

where the second equality applies to the constant density case. Clearly, with A = the mean velocity 
cannot become zero. Let us contrast it to the case with A > 0. One obtains 



2M 



I- A 



vac 



(70) 



As can been seen from (|7Uf) the mean velocity in the ellipsoidal configurations is decreased because of 
the K- external force. A drastic effect of the cosmological constant could be reached for the geometrical 
factor A approaching the critical value -4crit = p/pvac which is possible for very flat objects. In the 
extreme the mean velocity can go to zero ^ Together with inequality H16|l . the result in Eq. H70|l tells 
us that the temperature , T oc {v^) of the objects is very small if the objects is above the limit to reach 
equilibrium (the square bracket in Eq. (|70|) is then very small). This is a qualitative conclusion based 

^In the case of angular velocity, i.e Eq I54i this is somewhat different since f2o — > as e — > 1 independently of the 
density. 
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Figure 5: Behavior of Gmax as a function of the ratio pvac/pc for oblate (sohd hne) and prolate (dashed line). In the left 
hand side we show the complete range, In the right hand side the same plot for e — ► 1 and p -C Pvac 

in the presence of non-zero cosmological constant. For ellipsoidal configurations, this defines also a 
maximum value for the eccentricity emax given p or vice- versa, i.e, a minimum value for the density paun 
given e, through the relation pmin = -4(emax)Pvac- The bahavior of Cmax is shown in Fig. 5 as a function 
of Pvac I P- Galactic clusters with Pvac/p ~ 0.1 may have a vanishing mean velocity for e^ax ^ 1 hi the 
oblate case and Cmax ~ 0.92 in the prolate case. If the density is smaller, a vanishing mean velocity can 
be reached for non so flat objects. Spherical configurations have (u^) for p w 2pvac- 
Since the mean velocity is proportional to the temperature, p = M/V , Eq. H7U|I represents also a 
mass-temperature relation. Hence the results for the mean velocity squared are also applicable to the 
temperature of the configuration. For instance, we can write Eq. H69|l for a cosmological structure, 
say a galactic cluster, by writing its density as resulting from a perturbation 5p from the background 
density of the universe Ph{t) as p^c = Pb(i)(l + ^(O)j where 5{t) = 5p{t)/ ph{t)- Equation iTTUIl allows us 
to determine the temperature of the cluster at a given cosmic time t as 



He) 




1 - Aie) 



a + mr 



(71) 



arcsm e 



Prolate, 
Oblate, 



where a{t) is the scale factor and p is the mass of average components of the cluster. Equation H71|l 
assumes a flat universe i^mattcr + f^vac ~ 1- This result is a generalization of a result derived in |29| 
valid for spherical geometry. In that case one recovers the typical mass temperature relation T cx M^/^ 
mantaining p constant. Although this has the same dependence as in equation (|30|) . the meaning is 
different since (|71|l for A = 2 computes the temperature of a certain galactic cluster at some redshift 
given its mass while equation (|30|) is associated to a reversible process where a configuration passes from 
T = to some final T^, through states of virial equilibrium mantaining a constant mass. 



5 Small oscillations in the Newton-Hooke spacetime 

The stability condition of Newtonian configurations against oscillations can be also derived from the 
second order virial equation by expanding the periodic Lagrangian perturbations ^(r, t) = ^(r)e"^* with 
oscillation frequency tu. For simplicity, we will consider a rotational configuration without internal 
motions. By assuming adiabatic perturbations, the variational form of the second order virial equation 
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® is written as (see |77lliniET] l 

- + y TTpvac) P (r.C, + e.^;.) d^' = p [ed.^l + (1 - r)<5,fee^a,$N] d'V (72) 

+ ^-^^fc ^(5 - 3r)p [(O • 0(f^ • r) - n^{v ■ 0] dV, 

where F = (SlnP/Slnp)^, is the adiabatic index governing the adiabatic perturbations. For polytropic 
e.o.s we have 7 = F, of course. The Lagrangian displacement is constrained with the boundary conditions 
z), ^ = at r = and ii) ^ must be finite at the surface. The boundary conditions are required to solve 
the equation of motion for ^ resulting from perturbing Eulcr's equation (the Sturm-Liouvillc eigenvalue 
equation) . 

Let us consider radial oscillations which are adequate for spherical systems and for small deviation from 
spherical symmetry. By assuming a general trial function for the Lagrangian displacement ^ = /(r)f 
satisfying the boundary conditions and taking the trace in Eq. H72|l . we can solve for frequency of the 
oscillation about equilibrium 



2 ^ ^ 
A 



r - ^ _ 2 ( F - ; 1 B - 



(73) 



where we have assumed F — constant throughout the configuration and we have defined in analogy to 
Eq. (|16() the quantities 

In Eq. (|74(l g is a parameter with units of density, and 



W^ = - / p(r)/(r)^dV, R=^nl,J pir)f{r)r-'[r^-S,,x,Xj]d\, I=j p(r)/(r)rdV, 

(75) 



together with B = R/\W\. The critical adiabatic index is written as 



Fcrit = ^ (1 - 2B) ' 



2 Q 2 



(76) 



such that for F < Fdit instability sets in and the system becomes unstable while perturbed. From 
Eq. (|73|) . we see that through the inclusion of the cosmological constant we (the system) are forced to 
choose a bigger adiabatic index. In the simplest case, when /(r) = r, we have W W, I ^ I, B ^ (3. 
Furthermore, for p ^constant, we get A A and g = p. The stability condition then becomes 



F> ^(1-2/3)-^ 



1 + \aP^ h 
1 p L 



(77) 



The spherical symmetry which we assumed compels us to write the ratio /? for low eccentricities or to 
parametrize it in terms of the total angular momentum of the configuration. In the first approximation, 
at low eccentricities the ratio (3 given in (|58|) takes the form 



^.M^ P 



vac 



3 V5 p 



^e' + 0(e4). (78) 



This equation is useful if we want to calculate a small, but peculiar effect. Insisting that F is very close 
to 4/3 we can convert the stability condition F > Fcrit h^to a condition on eccentricity, namely 

\ 1/2 / A \ 1/2 

- /^vac o \ I I Pvac 



e > EE 3.35 (^1 + 5^ - 3ej + ) - (79) 

where e = | — F measures a small departure from the critical value 4/3. For the specific case e = 0, we 
conclude that the eccentricity must be such that e„ii,i < e ^ 1 with 

^ 1/2 

Pvac 



= 3.8 ( ^ ) , (80) 
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Figure 6: Critical adiabatic index as a function of /3 for a homogeneous spherical configuration for different values of 

C = 2pvac/P 

which is clearly vahd only for large densities, say for p > lO'^Pvac- This imposes a range for the 
eccentricity Ccrit < e ^ 1 for stability under radial perturbations. 

For completely spherical configurations with constant density, we can write /3 = fl^^^/Airp. Figure [S] 
shows the behavior of the critical adiabatic index as a function of f3 for different values of C = 2pvac//0- 
The largest deviation is for low densities as expected. 
For non rotating configurations, the stability criteria is simply given by 



If we work with a constant density profile, the choice of another trial function satisfying the boundary 
conditions mentioned before will not play a role in the above expressions since A A — 2 for any f{r). 

Related results regarding the adiabatic index in cosmologies with non-zero A have been obtained in the 
relativistic frame work in f j'2 j . 

6 Conclusions 

In the present work we investigated in detail the astrophysical relevance of the cosmological constant for 
equilibrium configurations. Using the tensor and scalar virial equations and the Lane-Emden equation 
we could show that many astrophysical facets get modified by A. The second aspect concerns the fact 
that A introduces new relevant scales (these scales would be zero if A = 0) like the maximal virial 
volume defined by the maximal virial radius H24|l and the maximal extension of bound orbits given in 



It is often assumed that superclusters with densities ~ pcrit are not in equilibrium. With the inequality 
(|16|) we have a precise tool to quantify this statement. Indeed, the pancake structure of the superclusters 
lead us to the conclusionthat they are even far away from the equilibrium state due to the factor A 
which grows with the object's flatness. On the other hand, relatively low density objects can still reach 
equilibrium (even ii A = 2 for the spherical case, the density of the object has to be only twice as large 
as Pvac, i-e, 1.4pcrit)- There is nothing which could, in pronciple, prevent a relatively low density object 
to be in equilibrium This is not only a result obtained from the virial equations with A, but follows also 
from hydrostatic equilibrium with the inclusion of A |12| . It is then natural to put forward the question 
how the shape of such objects affect their properties as compared to the case without A. For instance, 
if p/pvac — 10 and e ^ 0.9, the angular velocity is reduced approximately 30% with respect to its value 
with A = (firot = 0.7 flo). Other effects discussed in the text are the mean velocity (see Fig. 5) and 




Stt p{r)r'^f{r)dr 
3 '^Cpiryfir)^dr 



(81) 



Eq. (EH). 
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the adiabatic index (see Fig. 6). To be specific let us takje an example of a very low density object, 
p = 2.566pcrit- This object is still virialized as long as its eccentricity is e < 0.5. Approaching e 0.5 the 
mean velocity (temperature) goes to zero until finally beyond 0.5 the objects ceases to be in equilibrium. 

To highlight some quantities which get affected let us mention the polytropic index n and the 
angular velocity. Low-density objects with p ~ 20pcrit an index n > 1.5 becomes unphysical as the the 
conglomeration of matter does not have a definite finite radius. Taken together with the results of section 
five on stability against small oscillations this implies that the adiabatic index 7 from the equation of 
state P (X p^ gets restricted for these densities in a rather strong way, namely by 1.33 < 7 < 1.66. 
For higher densities the effect grows with n. For instance, for p ~ lO^pcrit, n = 5 is an unacceptable 
solution due to the lack of a finite extension. For non-spherical configurations we could show an effect 
of A on angular velocities, the mean velocity of the component of the astrophysical object and on the 
critical mass, again for densities one and two orders of magnitude above the critical one. Wc think that 
the work can be generalized in various ways. For instance, to generalize the virial equations to be able 
to differentiate between different models of dark energy mentioned at the beginning. Another aspect is 
to look into the Lane-Emden equation for non-spherical geometries ^1 02] • Here we concentrated on 
the equilibrium condition of already virialized matter. The collapse of matter in an accelerated universe 
with dark energy has been considered in |44) . 

Appendix 

Here we briefly show the solution for the cubic equation (|19|) . Let us write that expression as 

y'^+Vy + q = Q, p=10?7ri>0, g = -3r,ri < 0, (82) 

corresponding to a positive cosmological constant. Associated with this expression one defines a dis- 
criminant and a quantity R defined as 



^ = ^p' + ig2>0, R^sgn{q)sj\\p\<Q. (83) 
Since p > 0, the roots of (|82|) are given in term of the auxiliary angle 4> defined as 

sinh<^=^, (84) 
so that the only real solution for H82(l is written as 

" (85) 



y = -2i?sinh ( ^ j = ^ sinh 



iarcsinh f ^) 
3 \2R^) 



The other two roots are complex numbers whose real part are negative. If p < and q > Q the roots 
depend whether D is bigger or smaller than zero. If < 0, the three roots arc real and negative, while 
D > Q yields one real negative root and two complex roots which real parts are given as 



\l ^|p|cosh 



-arccosh f — ^1 
3 \2R^) 



(86) 
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Abstract 

Using the scalar and tensor virial equations, the Lane-Emden equation expressing the hydrostatic 
equiUbrium and small oscillations around the equilibrium, we show how the cosmological constant 
A affects various astrophysical quantities important for large matter conglomeration in the universe. 
Among others we examine the effect of A on the polytropic equation of state for spherically symmetric 
objects and find non-negligible results in certain realistic cases. We calculate the angular velocity 
for non-spherical oblate configurations which demonstrates a clear effect of A on high eccentricity 
objects. We show that for oblate as well as prolate ellipsoids the cosmological constant influences 
the critical mass and the temperature of the astrophysical object. These and other results show 
that the effect of A is large for flat astrophysical bodies. 

keywords: Large-scale structure of Universe. Galaxies: clusters: general. Instabilities. 

1 Introduction 

It is by now an established fact that the universe accelerates faster than previously anticipated 
000. Hence some hitherto neglected ingredient (in general called Dark Energy) has to be responsible 
for this new phenomenon. To account for this phenomenon we can introduce new physics in terms of a 
scalar field Q or modify three expressions in Einstein's equations which are often equivalent to a specific 
model with a scalar field. The three possibilities to account for the new physics are: the Einstein tensor, 
the energy-momentum tensor (this is to say, the energy momentum tensor of a fluid gets modified by 
the inclusion of other components, 3; or the equation of state 0). The first possibility encompasses 
a positive cosmological constant A and higher order gravity with a more complicated Einstein-Hilbert 
action In the the present work we choose to work with the cosmological constant as the simplest 

explanation for the acceleration of the universe. We shall put forward the question if such a cosmological 
model has an influence on astrophysical structures. We shall use equilibria concepts like hydrostatic 
equilibrium and virial equations to see how relatively low density astrophysical matter of different shapes 
behaves in a fast expanding universe. Anticipating the results, we can say that indeed there are some 
interesting effects. 

Often it is assumed that A does not have any effect on astrophysical processes which take place at scales 
different from the cosmological ones. Indeed, looking at the scales set by A, this assumption seems to be 
justified at the first glance. The scales set by the cosmological constant are of truly cosmological order 
of magnitude ^3^3)- The density scale is set by A = STrGNPvac with pvac — (0.7 — 0.8)/9crit- The length 
scale, r\ = 1/\/A is of the order of the Hubble radius while the mass scale M\ = r\/GN reaches up to 
the value of the mass of the universe. These scales constitute the so-called coincidence problem, namely 
the question as to why we should live exactly at an epoch where the scales of the cosmological constant 
are also the scales of the universe. Neither was it so in the past nor will it be so in the future when the 
universe expands further. The only astrophysical structures which match these scales are superclusters 
whose densities are indeed of the order of magnitude of pcrit- Indeed, here we can almost be sure that 
the cosmological constant is of relevance ^3^3- However, probing into astrophysical consequences 
of the cosmological constant of other, smaller and denser structures like clusters of galaxies or even 

* E-mail: a-balaguC!)uniandes.edu.co 
^ E-mail: mnowakosC!)uniandes.edu.co 
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galaxies themselves, would look a hopeless undertaking unless we find circumstances where the efli"ect 
of A (which in the very principle is present) gets enhanced. This can indeed happen through various 
mechanisms. For instance, in a problem where ta combines with a much smaller length scale, say ro, 
the efli"ects can be sometimes expressed as rgr™. In consequence, the observable in which this expression 
enters gets afli"ected by A in a way which is important at much smaller scales than ta. A concrete 
example is the Schwarzschild-de Sitter metric where we find the parameter ta together with the much 
smaller length scale of the Schwarzschild radius r,,. These two conspire in the form (rgr^y^^ to define 
the largest extension of bound orbits as explained in the text and in We will discuss a very similar 
combination which emerges from the virial theorem defining the largest possible virialized structure 
with a given mass. Another possibility to enhance the effect of the cosmological constant is to consider 
non-spherical objects. It then often happens that the effect of A becomes {h /h)'' Pvac where are two 
different length scales of a fiattened object like a disk or an ellipsoid This indeed happens for many 
astrophysical quantities, among other the critical mass, the angular velocity and the temperature (mean 
velocity of the components of the large structure) which we will discuss in the present paper. Finally, 
we can vary a dimensional variable to see if this enhances the effect of the cosmological constant. As 
an example let us quote the polytropic index n in the equation of state entering also the Lane-Emden 
equation. It is known that with growing n (n > 5) the object described by this equation of state does 
not have a well-defined radius as the density goes only asymptotically to zero. We will show that this 
pattern of behaviour becomes more dominant with A ^ 0. 

Of course, all these effects become stronger the more diluted the mass conglomeration is. The superclus- 
ters are certainly the best candidates if we look for astrophysical effects of the cosmological constant. 
As a matter of fact, they do not seem to be virialized due to the extreme low density and their pancake 
structure ^3 where the effect of fiatness mentioned above becomes powerful ^3^3- For the next 
structure, the clusters of galaxies (or groups of galaxies) with densities between one and three orders of 
magnitude above the critical density we would need one of the enhancing factors discussed above 
to see an appreciable effect of A. This is possible in various ways as shown below. Galaxy clusters can 
have various forms, among others oblate and prolate ^3^^. And what is more, they can even rotate 
We will show an explicit effect of A on their angular velocity and temperature in case the angular 
velocity is zero. Since the effect of large eccentricity is larger for prolate than for oblate ellipsoids, 
it is comforting to know that clusters can assume a prolate shape. For low-density galaxies like the 
Low Surface Brightness (LSB) galaxies whose density is roughly four orders of magnitude above the 
critical density ^3^3; we still find some effects. For n = 5 and A = the solution of the Lane-Emden 
equation is very often used as a phenomenologically valid description of the density profile (called 
also Plummer's law This is still possible as p — as r — cxd. However, this property vanishes for 
low-density galaxies and the n = 5 case not only does not have a well-defined radius, its solution does 
not vanish asymptotically which thus rendering it unphysical. 

The paper is organized as follows. In the second section we will briefiy review the general form of 
virial theorem including pressure, magnetic fields and, of course, the cosmological constant. Here we 
will also discuss some general results regarding A. In the third section we will specialize on spherical 
configurations. We will show how A sets the scale of a maximal virial radius and compare it to a 
result from the Schwarzschild-de Sitter metric. We will also solve the Lane-Emden equation numerically 
and analytically (for the polytropic index n = 2). In the fourth section we will discuss non-spherical 
configurations. First, we will show how A affects the angular velocity of spheroids. In addition, we will 
discuss the effects of A for the critical mass, mean velocity and mean rotational velocity too. The fifth 
section is devoted to small oscillations around equilibrium. 

2 Local dynamics with cosmological constant 
2.1 Newton-Hooke spacetime 

The cosmological constant enters the equations of Newtonian limit as a consequence of its appearance 
in the Einstein field equations. It is through this weak field limit approximation that A enters also in 
the equations describing the structure of astrophysical configurations. It is interesting to note that all 
variables to be found in the virial equations, are also present in the Poisson equation of the Newtonian 
limit. However, this is not always the reason why these terms enter the virial equations, at least in the 
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first order. The Poisson equation for a self-gravitating system modeled as an ideal fluid is written as 



V2<D = 4^GAr(p + 3^ + 2^) -A. (1) 



where P is the pressure and Ueva is the electromagnetic energy density. The solution of Q at the zeroth 
order of v / c (from now on we set c = 1) is written as 

<D(r) = -Gn I d^' - -k\v? 



v 



|r — r'l 6 



where the dots stand for the correction terms that appear because the boundary conditions are now 
set at a finite distance These terms can be usually neglected. The cosmological constant A > 

contributes to the expansion of the universe. This fact remains partly valid in the Newtonian limit 
where A gives us an external force. This defines the so-called (non-relativistic) Newton-Hooke spacetime 



2.2 The A-virial theorem 

The second order tensor virial equation can be derived in different ways: from a statistical point of view 
through the collisionless Boltzmann equation, from a variational principle or by direct differentiation of 
the moment of inertia tensor 

p{r)rirkd^r. (3) 

In the following we use the statistical approach ^3 which also allows to derive higher order virial 
equations (for instance, the first order virial equation refers to the motion of the center of mass). In this 
context, from Boltzmann's equation one can derive the equation for momentum conservation (Euler's 
equation) written for a self gravitating system infiuenced by a magnetic field as 

+ M<D + \d,{B') + d, {V,, - B,B/) = 0, (4) 

where $ is the gravitational potential given by Q (which includes A) and 

Vij = p{{vi - {vi)){vj - (vj))) = SikP + TTik, Tr(7rife) = 0, (5) 

is the pressure tensor, P is the pressure and TTik its traceless part. Equations Q and ^ together with 
an equation of state P = P{p, s) [s is the entropy) complete the description of a self gravitating fiuid. 
By taking exterior products of with Euler's equation and integrating over the volume of the system 
one obtains the second order virial equation as 

- 2Tife - I Wf™ I -I- -TrGAT/Ovac life + Ilife where = / Vikd^r, (6) 



2 dt2 I »fc I ■ 3 



V 



where Tik is the kinetic energy tensor and iVVf^*^"] is a generalized potential energy tensor which contains 
the contribution from the gravitational potential energy tensor Wfj. and the contributions of magnetic 
field through 

|Wri = |Wr,|(l-A(,,)). (7) 
The other quantities are defined as follows 

= \j p{v.){vk)d\, Wl = -Gn f p{r)rAmdh; A(,,) = (8) 

^ Jv Jv KVifel 

together with 



TikiB) = SikB - 2Bik - rk 

JdV 



1 2 

-SijB - BiBj 



dSj , Bik = \ Bi Bk d^r, B = lY{Bik ) . (9) 
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A very useful version of the virial equation can be derived by assuming an isotropic pressure tensor and 
taking the trace in This way we get the scalar A- virial equation 

i0 = 2/C - |W8-| + ^t:GnP..cI, (10) 
where the total kinetic energy is written as 

K = - ( p{v^)A^r = T with n = / PA^r (11) 



The equilibrium condition is reached for 1 = 0. This gives us the general A- virial theorem 

2Tife-|wf"| + ^^Pvac2:ife + nij = o, 2/c- |W8^"| + ^^Pvac2: = o. (12) 

For rotating configurations with constant angular velocity, the kinetic term is modified as in the standard 
way as 

Tik ^ Tik +TZik, T^ik = i^Tot^ik — ^Totilkj^)"^) , R = Ti(TZij), (13) 

with TZij the rotational kinetic energy tensor and Tik is referred to motions observed from the rotating 
reference frame. The A-virial theorem has been used in different contexts in In the 

present work we will extend these studies. 



2.3 General consequences 

The tensor virial equation is widely used in many astrophysical applications. The inclusion of A provides 
a new way to study effects of the cosmological constant (parameters , in general) on astrophysical 
objects. The outcome depends essentially on two factors: the geometry of the configuration and the 
density profile. We will explore the spherical geometry for both constant and varying density profiles 
and study some effects for non spherical geometry with constant density. The consequences that can 
be derived from the A-virial theorem can be classified in two categories. The first one puts an upper 
bound on the cosmological constant or alternatively a lower bound on density of objects in gravitational 
equilibrium. Provided these bounds are satisfied, we can also study in the second step the effects of A 
on other properties of the astrophysical configurations like rotation, small oscillations etc. 
The first simple consequence of the virial equation emerges if we require the system to satisfy ^3 • The 
fact that /C > implies an upper bound on the vacuum energy density 

All systems in equilibrium have to satisfy . Note that the right hand side of this expression is a 
function of both the density and the geometry of the system. Hence, we must expect different bounds 
for different geometries and density profiles. For instance, if we assume a constant density and B = 0, 
we can define and I through 

|w^| = iGVl>v^l, i = pi, (15) 

such that the bound written in becomes 

,>^,_, with A^'-f(^-^^. (16) 

The factor A which is only a function of the geometry (if we neglect the contribution of magnetic 
fields), will appear in many places in the paper. Its relevance lies in the fact that it enhances the effect 
of the cosmological constant for geometries far from spherical symmetry when A is large. A useful 
generalization can be done for situations in which we use the tensor form of the virial equation, namely 

_ IGtt / lij \ 

An = — ^ — , p = constant. 17 

' 3 [\Wfp\)' ^ ' 
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Finally, a curious equation can be derived by eliminating A from the tensor virial equations: 

27^i-|nf"l + n., 

2r - iw8™i + n ~ T ' ^ ' 

Although A does not enter this equation, ^3 is only valid if the denominator is non-zero as is the case 
with A 7^ 0. In Sect. 4 we will use this equality to infer a relation between the geometry and rotational 
velocity of an ellipsoid. Having discussed the general form of the virial theorem, we will discuss now the 
effects of A and set B = 0. 



3 Spherical configurations 

The tensor virial equation is trivially satisfied for spherically symmetric configuration without a magnetic 
field, since W/^ = i^ifeW^ and lik = 5ikl. Therefore, in this section we use only the scalar form of 

3.1 Constant density 

Explicit expressions can be derived in the spherical case with constant density, with |>V^| = f '^'^^ 
and X = |Mi?^, so that ■4sphericai = 2. In this case the ratio Pvac/p does not get enhanced much by the 
geometrical factor A. 

Worth mentioning is the result from general relativity. There the upper bound for the cosmological 
constant comes out as A < AttGnP where p is the mean density defined hy p = M/V. This 

bound is derived not only from Newtonian astrophysics, but also from a general relativistic context via 
the Tolmann-Oppenheimer-Volkoff equation ^3 for hydrostatic equilibrium of compact objects ^3^3- 
Another relevant effect of the cosmological constant is the existence of a maximal virial radius of a 
spherical configuration which can be calculated from the A- virial equation. Using the expressions for 
jW^I and X given before, equation yields as a cubic equation for the virial radius i?vir 

Rl^^ + (lOryri) R,,, - 3r,rl = 0. (19) 

Here we introduced the dimensionless temperature rj as 

_ /C 3^5 

V = — = T, 20 

where p is the mass of the average member of the configuration, ks is the Boltzmann constant, T is the 
temperature and is defined by 

= GnM, (21) 

The length scale r\ is set by the the cosmological constant as 



TA = = 2.4 X 10^ /ifo^ I X iQi" ly. (22) 

hro 0.7 is dimensionless Hubble parameter and Slyac = Pvac/Pcrit ^ 0.7 is the density parameter 
at the present time. The positive real root of equation is given by 

i?™(??) = C^(??)i?vir(0), (23) 

where i?vir(0), is radius for the configuration at ?y = is given by 

i?vir(0) = (3a;)'/^ TA = (3r,rA)l/^ (24) 
and the dimensionless parameter x is defined as 

r f M \ 

^=^ = 1-94 X 10-23 j^—j hro^li', « 1. (25) 
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The radius i?vir(0) is the largest radius that a spherical homogeneous cloud may have in virial equilibrium 
(i.e, satisfying ^3). The function zu(ri) can be obtained from the solution of the cubic equation and 
reads (clarified in the appendix A) 



w{ri) = 2.53a;-^/^?y^/2sinh 



iarcsinh (o.^Axr] ^1'^ 



(26) 



Figure 1 shows the behavior of n7(?y) for different values oix. We see that the increase of the temperature 
implies a decrease of the effects of A which can be easily checked if we solve i?vir from the virial theorem 
with A = and compare it to the approximation ?y — cxd in : 

i?vir(A ^ 0) = i?vir(?? ^ OO) = i?vir(??*) = ^ — • (27) 

We can consider as a radius-temperature relation for a fixed mass applied on astrophysical structures 
in a single state of equilibrium in the presence of A. That is, given x and rj we calculate the radius. But 
we can adopt another point of view for this relation. Imagine a spherical configuration characterized by 
a constant mass M . In analogy to a thermodynamical reversible process, the configuration may pass 
from one state of virial equilibrium to another following the curve w — rj, that is, satisfying the condition 
1 = 0. Clearly, there must be some final temperature rj^ when this process ends since the temperature 
cannot increase indefinitely. But of course since the virial equations are not dynamical we cannot know 
which stage is the final one. If we assume that the effects of A are negligible when rj = rj^, then using 
Eqs ESI and El 

we get 

n ono -1 2/3 -Rvir('?*) /„on 

?y* = 0.208ro^ x' , ^* = -5 — 77^- (28) 

-ftvir(.Uj 

This is an equation for the temperature 

M-2/3t^ = 0.138 (^^^ (29) 

For a hydrogen cloud (jj = mproton), we then write the mass-temperature relation using Eqs.^^ and 
^3 as 

/ (iJ \ 2/3 

n = 8.60 X W-'w;' ^—j h^'^Qlil K. (30) 

Note that this expression maintains the same dependence of the standard mass-temperature relation 
derived from the virial theorem, i.e, T cx M^/^ (see ^3 or equation ^3 of this pap er). However the 
meaning of is different from that of typical mass temperature relations since is associated 

to the temperature that a system acquires in the final stage after after going through some reversible 
processes which took the system through successive states of virial equilibrium with constant mass from 
a radius i?vir(0) to a radius i?vir(??*) or vice versa. On the other hand, the mass-temperature relation 
like Eq. ^3 of this paper relates the temperature of any configuration in equilibrium with its observed 
mass at constant density. In that context one considers only one equilibrium state and the cosmological 
constant enters just as corrections. 

As a final remark on Eq. we discuss a result which formally coincides with the virial radius i?vir(0) 

derived from the Schwarszchild de-Sitter spacetime i fi , but whose physical meaning is quite different. 
The Schwarzschild-de Sitter metric takes the form 

ds2 = _^_^_Z!_^dt2+ fl-^-^^ dr2 + r W + r2 sm2 (31) 
\ r 3r^ / \ r 3r^ / 

Now in contrast to we have 

Rs = Gn//, (32) 

with 1^ the mass of the object giving rise to the Schwarzschild-de Sitter metric (in Eq. M is the mass 
of the total conglomeration whereas here we consider /j, as the mass of its average member). Choosing 
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Figure 1: Ratio between Rviv{v) SLnd i?vir(0) for different values of :r as a function of = (3//i)T, where /i is the mass 
of the main average components of the system. 



the afEne parameter as the proper time r the equation of motion of a test-body can be cast in a form 
similar to the corresponding equations from non-relativistic classical mechanics. 

+ [/efj = i (5^ - 1) = C = constant, (33) 

where 5 is a conserved quantity and (Jeff is the effective potential, defined by 

d/ R 1 

5=(l + 2[/eff(r))— , UMr) = -—-^-- (34) 
dr r 6 

For simplicity we are have chosen here the angular momentum L to be zero. With L zero or not, (Jeff 
displays a local maximum below zero due to A 7^ forming a potential barrier. This is to say, the 
standard local minimum where we find all the bounded orbits is now followed by a local maximum after 
which (Jeff goes to — cxd. With A = this function approaches zero asymptotically. One is immediately 
tempted to say that this barrier will occur at cosmological distances. This is not the case and one 
calculates 

r^,, = imsri)''' = 9.5 x IQ-^ ( j^Y' ("^V Mpc. (35) 

In other words, the combination of the large scale r\ with the small scale Rg gives us a distance of 
astrophysical relevance, namely r^nax- Its relevance lies in te fact that beyond r^nax there are no bound 
orbits. Indeed, with /j, the solar mass, r^nax is of the order of a globular cluster extension (70 pc); with 
/J, as the mass of globular cluster, r^nax comes out to be of the order of the size of a galaxy (10 kpc), 
and finally taking /j, to be the mass of a galaxy, r^nax gives the right length scale of a galaxy cluster (1 
Mpc). Certainly, the value of the extension of a large astrophysical body is the result of a multi-body 
interaction. But with the actual values of r^nax, it appears as if the length scale (we emphasize that 
we are concerned hare about scales and not precise numbers) of an astrophysical conglomeration is 
approximately r^nax, which apparently means that this scale does not change drastrically when going 
from a two body problem to a multi-body calculation. This makes sense if the object under consideration 
is not too dense. We can now say that whereas M in i?vir(0) (via r^) is the mass of the object, /j, in Eq. 
ESI is the mass of its members. Clearly, we have i?vir(0) ^ fmax, but both scales are of astrophysical 
order of magnitude. A result related to derived in the framework of general relativity can be found 
in 



1 /dr 

2 Vd7 



3.2 Non-constant density 

The examination of configurations with non-constant densities can be done in two directions. Knowing 
the density profile p{r), we can set up the virial equation and evaluate the equilibrium conditions from 
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the inequality In this picture, the effects of A are included in the solution for the potential $ as 

in Eq. Q and the resulting term acts like an external force, as mentioned before. 

The second option is to combine the Eqs. Q, ^ and an equation of state (e.o.s) for which we can take 

1+1 



a polytropic form P 

e V de 



Kp 



r 



where pc is the central density 
Jeans length defined as 



. This way we obtain the Lane-Emden equation with A ^3 
Cc Cc = 2 (^^) , (36) 
= a^, p{r) = Pc'tp" {£.) with tp{0) = 1, ^"'(0) = ^nd a is the associated 



lK{n + 1) 



Airpc 



(37) 



It is important to notice that in this picture the expected effects of A are to be found in the behavior 
of the density profile since now Eq. ^3 implies that its solution is also function of the parameter Cc- 
Some effects of A are contained in the total mass and the radius of the configuration which is reached 
when ip(^i ) = 0. Hence implies 
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1/2 



R = a^i 



rA 



K(n + l)pc 



1/2 



(38) 



Note that the radius of the configuration is now proportional to r\. This is due to the fact that A 
sets a scale for length. However, this does not mean that R will be always of the order r\ as A is also 
contained in the expression in the square brackets in Eq. We will show this below in a concrete 

example. Since A is a new constant scale the Lane-Emden equation loses some of its scaling properties 
as explained in The mass of the configuration can be determined as usual with. 



M(0 



Aira^Pc 



de 



-Tra^^Vc 



Cc 



(39) 



IdC\ 

where we used Eq. ^3 for the second equality. The total mass is then obtained by evaluating the last 
expression at ^ = ^i. As expected, the mass increases because the Newtonian gravity has to be stronger 
in order for the configuration to be in equilibrium with A ^ 0. Figure 2 shows the numerical solutions 
for n = 1 to n = 5. We expect that the radius of the configuration is increased by the contribution 
of Cc and find it confirmed in the figures. However, not always is the radius of the configuration well 
defined, even if n < 5. For sizeable values of Cc (black line) we cannot find physical solutions of Eq. 
^3 as the function ip acquires a positive slope. One might be tempted to claim that the radius of the 
configuration could be defined in these situations as the position where ip has its first minimum, but as 
can be seen for n = 3 such a radius would be smaller than the radius with Cc — which contradicts the 
behaviour shown for the other solutions where ^(Cc 7^ 0) > ^(Cc = 0). As already mentioned above this 
is the correct hierarchy between the radii because large Cc gives rise to a large external force pulling at 
the matter. The numerical solutions show that for relatively large values of Cc, only n = I has a well 
defined radius. In this case the effect of A is a 13% increase of the matter extension as compared to 
= 0. As we increase the polytropic index, Cc ^ 10~^ leads to non-physical solutions while the effect 
with bigger values of Cc becomes visible only for n = 3. For instance, Cc ^ 10~^ results in a radius 
which is 17% bigger than the corresponding value with Cc — 0. The combination n = 4 and Cc ^ 10~^ 
also leads to a non-physical solution, whereas the radius of the case Cc ^ 10"'* displays a difference of 
13% as compared to Cc 0. Finally, for n = 5, the only physical solutions are obtained for the lowest 
values of Cc where ip' < This case is particularly interesting as with A = it is often used a a viable 
phenomenological parametrization of densities ^3^3- The solution has an asymptotic behaviour as 
r~" which has been also found in LSB galaxies de Blok at al. 2004). With A 7^ the n = 5 seems less 
appealing as the matter is diluted. For all values of n, the difference between Cc = 10~® and 10~^ is 
negligible. 



Analytical solutions of ^3 can be found for n = 0, 1 and n = 5 if A = 0. For n — 
equation of state reduces to the equation of state of the isothermal sphere P = Kp. 
A 7^ 0, we can write the analytical solution in the case n = 1 as 



► cxD, the polytropic 
As an example, for 



fi 



Cc)^+Cc 



(40) 
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Figure 2: Effects of A on the behaviour of the density of a polytropic configuration for different ratios C^c and different 
polytropic indices. The radius of the configuration is not always defined, even for n < 5. For higher values of Pvac , only 
the n = 1 case has a definite radius for these values of C^c- For other cases, the configuration is defined only for small C^c- 



The radius is R = a^i , where is the solution of the transcendental equation 

iiCc = -(1 - Cc) sin^i. 
In the first order of Cc one finds 



R. 



TTK (1 + Cc) 



(41) 



(42) 



Equation ^3 ^Iso implies that there exists some Ccrit such that for Cc > Ccrit, we cannot find a real 
solution for . Approximately this gives 



Pc > lO.Spv 



(43) 



which, provided the overall density is not too big, is better that p > 2/9vac which is a result from the 
general inequality for p = const and spherical symmetry. Finally, we can calculate the contribution 
of A to the total energy of the object. Generalizing the results found in ^3 we obtain 



E=-(3-n)W'', 



-G 



N- 



R 



(44) 



from which we infer that the correction is very small in this case. 

At the end of this section we would like to summarize the findings from Fig. 2. In table 1 we write the 
ratio ^i(A = 0)/^i(A 7^ 0) for the same ratios Cc and polytropic index as in Fig. 2. The horizontal line 
represents a non-defined radius. The symbol cxd indicates that the radius is defined only asymptotically 
in case of A = 
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Cc 


n = I 


n = 3/2 


n = 2 


n = 3 


n = A n = 5 


0.1 


~ 0.88 










0.005 


~ 1 


~ 1 


~ 0.98 


~ 0.86 




2 X 10-4 


~ 1 


~ 1 


~ 1 


~ 1 


-0.88 


1.2 X 10-5 


~ 1 


~ 1 


~ 1 


~ 1 


~ 1 


6.1 X 10-^ 


~ 1 


~ 1 


~ 1 


~ 1 


~ 1 CXD 



Table 1: Values of the fraction ^i(A = 0)/^i(A 7^ 0) for different values of the ratio C^c and the polytropic index. The 
horizontal lines represents the non well defined radius. 



4 Nonspherical configurations 

4.1 Rotating configurations 

As emphasized already before, the effect of A can get enhanced for non-spherical objects. This happens 
when the vacuum energy gets multiplied by a ratio of two length scales li and I2 and we end up 
with expressions like Pvac{h/h)" ■ For instance, for ellipsoidal configurations the geometrical parameter 
A entering among others the inequality can be calculated from its definition ^3 with Xij = 

g7raia2a3(5jjaf and Wik given in We have O 



4 



A 



4 / 3 - \ e ai>a3 8 f ai 



3 yarcsine/ 2Vl — ^tt yas 



4 e(3 - 2e 



2\ 



pro 



3 (1 - e2)3/2 



In 



1 + e 
1 - e 



-1 o / 



3 \ai 



2a3A 



In 

ai J 



(45) 



where the eccentricity"'^ e is = 1 — flg/af for the oblate and = 1 — af/og for the prolate case. The 
behaviour of ^obi and ■4pro is shown in figure ^3 It is clear that for very fiat astrophysical objects 
we can gain in this way several orders of magnitude of enhancement of the effect of A if ^ is a factor 
of attached to pvac. Needless to say that we often encounter in the universe fiat objects like elliptical 
galaxies, spiral disk galaxies, clusters of galaxies of different forms and finally superclusters which can 
have the forms of pancakes. Of course, the more dilute the system is, the bigger the effect of A. We can 
expect sizeable effects for clusters and superclusters, even for very fiat galaxies. Regarding the latter, 
low density galaxies like the nearly invisible galaxies are among other the best candidates. 
A convenient way to model almost all fiat shaped objects is to consider ellipsoids which in the limit of 
fiattened spheroids can be considered as disks. There are three different kinds of elliptical configurations, 
characterized by three semi-axes ai = a, 02 = h and 03 = c: oblate , with a = h < c, prolate with 
a = h y c and triaxial systems with a > h > c. Here the tensor virial equation provides a tool to 
determine which of these geometries are compatible with the virial equilibrium. Considering the case 
A = and p =constant, or spheroids with confocal density distribution whose isodensity surfaces are 
similar concentric ellipsoids ^^^3 i-c 



3 2 

p = p{m'), m' = alJ2^, (46) 

1 

the oblate ellipsoid (MacLaurin) emerges as a solution of the virial equations with a bifurcation point to 
a triaxial Jacobi ellipsoid . We can view the virial equation without A as a homogeneous equation. 
Switching on A 7^ this becomes an inhomogeneous equation whose right hand side is proportional to 
Pvac • It is therefore a priori not clear if in the case of A 7^ we can draw the same conclusions as with 
A = 0. 

Let us assume a configuration which is rotating with constant angular velocity around the z axis. 
Neglecting the internal motions, the A-tensor virial Eq. Q for such a configuration is 

firlt i^ik - Si.Izk) -\yvll\ + ^ttGnP^A = -Sikll- (47) 



^ Once the density is given, the inequality becomes in this way a defining equality for emax such that e < 
order to maintain equilibrium. 
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Figure 3: Function g{e) and geometrical factor ^(e) for prolate and oblate ellipsoids. These functions have their largest 
values for flat oblates and large prolates. 



The equations with diagonal elements yield the following identities 



8 



\Wy 



This set of equations can be resolved for the angular velocity: 



7 

vac-^yy 



Pvac-^zz ■ 



8 f Izz 
-TrLrNPvac I 1 " J 



(48) 



(49) 



The same expression holds if we make the replacement Wxx '^yy ^nd I^x lyy on the right hand 



side of 



Si 



We can also eliminate 51;:^^ from these expressions to arrive at the condition 



\-W^x\-\-Wl\ 



L|>VN|-|WN| 



87r 



GnPv 



1-Si 



|WN|-|WN 



(50) 



with Si = lyy/Ixx- This expression determines the possible ellipsoidal configurations. One solution is 
obvious: Si = 1 together with |W^j,| = |VV^|. This happens if the density is constant or for a confocal 
ellipsoid with the density given in Hence this is not different from the A = case. However, the 

bifurcation point to a triaxial system will get affected by A Furthermore, the explicit calculation 

of this bifurcation point will now depend also on the details of the density profile even if we take the 
latter to be as in Eq. ESI- This is a direct consequence of a 'pre-existing' density scale pvac- For a 
homogeneous configuration we have 



Sii^-^ 

Tlx 



SttGnP 

Vx 



1 



— 1 A n 




^XX ^X^Z 



1 



rix 



(51) 



-%/ \ P J \^-% J 

where Si = Xyy/X^x, = l^i^ 1/1^/7 1 ^nd An are defined in Eq. For an oblate configuration they 

can be calculated in terms of the eccentricity e 



^1 — a^/c^ to give j^j 



Axx Ay 



3yr 



Az 



8 



37r 



Furthermore we have 

ry, = ryy = 2(1-6^)1/2 



arcsm e 



(52) 



(53) 
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Using these expressions, we write the constant angular velocity for an oblate ellipsoid from the first line 
of equation Eq. as 



Pvac 1 / X 

— ]9(e) 



where SIq corresponding to the angular velocity when A = is given by the Maclaurin formula 



"o 



8 



(l-e^) 



2x1/2 



-(3 — 2e^) arcsin e ^(l ~ 



(54) 



(55) 



The function g(e) defined through Eq. can also be calculated with explicit dependence on the 

eccentricity as 



f \ — -^xx Vx-^zz 4 g r 2n1/2/o o 2n • o /l 2\ 

g[e) = = -e (1 — e ) ' (3 — le ) arcsm e — 3e(l — e ) 

\ — Tjx 3 L 



(56) 



As is evident from the above equations, A has a twofold effect on the angular velocity. Firstly, it reduces 
the angular velocity with respect to the value SIq especially at the local maximum (see Figure 3). This 
is not a small effect and can affect even galaxies. Secondly, we see from Eq. that SIq — t- for e — 1. 
On the other hand 



p 



-9{e) 



Pvac 327r / 03 

P 9 v«i 



(57) 



approaching 1 for a very fiat oblate configuration and not too dense matter. Therefore, beyond the 
local maximum in Slpot the cosmological constant causes a steeper fall of Slpot towards 0. Another 
relevant interesting quantity which can be calculated in this context is the ratio of the rotational over 
the gravitational energy contributions to the scalar virial equations, i.e., 7?./|W^|. In accordance with 
Eq. ^3 the latter can be written as 



f3 = TZ/\W 



fi = fio 



p 



l3o 



3 
2^ 



e(l-e2)i/2 



(58) 



The effects on /3 are therefore similar to the the ones encountered in Slpot (see also Figures 3). Finally, 
on account of _P > we can infer from the virial equations with p =const the following inequality 



which together with results in an inequality for the density of an ellipsoidal configuration: 

1 fA-4:l3og' 



p > Apv 



A 



l-2/3o 



(59) 



(60) 



For e — 0, we have A — ?> 2 A, while for e 1, A ^ A, and A > A. Therefore the above inequality is 
slightly weaker than the bound given in . Nevertheless it is useful since it it derived directly from a 
different starting point [P > 0) than which is based on /C > 0. 



4.1.1 Mean mass- weighted rotational velocity 

Not always the deviation from spherical symmetry guarantees that the effect of A on observables is 
sizeable. This depends on the context and also which scales we compare. If we compare pvac to densities, 
the quantities A and An defined in Eq. are for fiattened objects large enough to enhance the effect 
of A. If r\ is combined with the Schwarzschild radius r,, to give r^^x as in Eq. the result is still of 

astrophysical relevance. But if we had to compare one of the axes of an ellipsoid to r\ (i.e. ai/r\) the 
effect would be negligible unless the extension a,- is bigger than Mpc (clusters and superclusters) and 
the small ratio ai/r\ is comparable to another small quantity of the same order of magnitude entering 
the equation under consideration. This happens for instance if we generalize a result (discussed ^3) on 
a mass-weighted mean-square rotation speed vq of an ellipsoidal object to include A. Assume that due 
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to symmetry properties of the object the only relevant components of the tensors in the tensor virial 
equations are xx and zz. We then obtain 



If the only motion is a rotation about the z axis we have T^^ = and we can solve for T^x as 



2Tx 



1 



X X ■J* ■J* I y\.^^'r» ■J* 
3 



(61) 



(62) 



Using Eqs. Q, Q and the quantities entering our equation can be parametrized in the following 



way 



27;. = pivDd^'r = l^Mvl, H.. = Mai, ^zz = (1 - <5o)n... 



(63) 



where Vq is the mass weighted mean angular velocity, ctq is the mass-weighted mean-square random 
velocity in the x direction and measures the anisotropy in Yin. If is of the order of one, it suffices 
to compare ctq with AI../3M. If both are of the same order of magnitude, the effect of A is non- 
negligible. Since vq and ctq are of the same order of magnitude and vq is non-relativistic, we can assume 
that (To < 10~^. The quantity A/../3M can be estimated to be {ai/r\)^. Hence if a,- ~ IMpc, vq (ctq) 
has to be truly non-relativistic and of the the order of 10~® to gain an appreciable effect of A. This 
improves if a,- is one order of magnitude bigger which is possible for large galaxy clusters. The velocities 
have to be then at most 10"'*. In these cases we have to keep A and while solving for one has 



Ivl 
2al 



(l-<5o) 



1 



1 Mr 



3 Mai 



Izz \W^x 



1 



(64) 



Note that if A cannot be neglected vljal is not only a function of the eccentricity for ellipsoids with the 
density give in but depends also on the details of the matter density as the latter does not cancel. 

Using almost the same set-up as above, we can use equation ^3 to establish a relation betweem^o, the 
mass M and the geometry of the object which we choose below to be oblate. After straightforward 
algebra we obtain 



2al 



3 Gf^M ^/Y^ 



3 arcsin e 



+ (i-M - 



1. 



(65) 



Since this relation is derived from ^3 which in turn is based on the assumption of A 7^ it is only valid 
for non-zero cosmological constant albeit the latter does not enter the expression. Note the enhancement 
factor (ai/as)^. 

4.2 Other effects of A for non-rotating configurations 

We can now derive other relevant quantities from the scalar A- virial theorem applied to homogeneous 
ellipsoidal configurations. In this section we will not consider rotating configurations, but systems with 
kinetic energy coming from internal motions. We will focus again on ellipsoidal, oblate and prolate, 
geometries. As in the preceding section the relevant quantity here is the function A written for both 
configurations in equation 

Critical mass: Consider the stability criteria for a homogeneous cloud with mass M and internal mean 
velocity (v) in mechanical equilibrium with the background with pressure P (see ^3). The system will 
collapse under it's own gravity (P < 0) if its mass is greater than a critical mass Mc- With cosmological 
constant, this critical mass is increased with respect to its value Mco when A = 0, which is expected, 
since now there is an external force associated to A that acts against Newtonian gravity and hence the 
collapse can be postponed. By using the scalar virial theorem ^3 and setting _P = as the criterion 
for the onset of instability, we can write for arbitrary geometry 



Mr = M, 
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1-A 



PcO 



M, 
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cO 
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(66) 
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Figure 4: Effects of A on the angular velocity O,^^^ and the ratio (3 for Pvac ~ 0.7pcrit 



This expression is valid for any geometry. However, we pointed out already that spherical symmetry 
implies A = 2 and the effect is suppressed. Some numerical values can be considered by writing 



SkBT/nip for a hydrogen cloud with T 500 K and radius R lOpc. One then has 



A 



Pvac 
PcO 



10" 



(67) 



which represents a very tiny correction to the critical mass Mco for a spherically symmetric object. On 
the other hand, for ellipsoidal configurations with the same temperature, we have 



SM = Aou '— 



Pvac 
PcO 



10" 



— y 



SM = A 



Pvac 
PcO 



10- 



' (kpc) 



(68) 



in the e — 1 approximation for oblate and prolate configurations, respectively. We have set c = 10a for 
the prolate case and Slyac = 0.7. For an oblate ellipsoid with a 50 kpc, one has SM Pd 0.15 while for 
the prolate ellipsoid with c Pd 50, SM Pd 10^. 

Mean velocity and Mass- Temperature relation. By virtue of the scalar A-virial theorem, we can also write 
down the mass-temperature relation for an astrophysical structure. Note first that using K, = ^M{v'^) 
in equation ^3 we have for A = the standard expression for the mean velocity 



/A=0 



\w 



Nl 



M 



M 



Nl 



where the second equality applies to the constant density case. Clearly, with A 
cannot become zero. Let us contrast it to the case with A > 0. One obtains 



(69) 



the mean velocity 



2M 



I- A 



Py 



(70) 



As can been seen from the mean velocity in the ellipsoidal configurations is decreased because of 
the A- external force. A drastic effect of the cosmological constant could be reached for the geometrical 
factor A approaching the critical value ■4crit = p/Pvac which is possible for very fiat objects. In the 
extreme the mean velocity can go to zero ^ Together with inequality the result in Eq. tells 

us that the temperature , T cx (v^) of the objects is very small if the objects is above the limit to reach 
equilibrium (the square bracket in Eq. is then very small). This is a qualitative conclusion based 



^In the case of angular velocity, i.e Eq this is somewhat different since 

density. 



as e 



1 independently of the 
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Figure 5: Behavior of emax as a function of the ratio Pvac/ Pc for oblate (sohd hne) and prolate (dashed line). In the left 
hand side we show the complete range, In the right hand side the same plot for e — >■ 1 and p <^ Pvac 

in the presence of non-zero cosmological constant. For ellipsoidal configurations, this defines also a 
maximum value for the eccentricity emax given p or vice- versa, i.e, a minimum value for the density pmin 
given e, through the relation pmin = -4(eniax)Pvac- The bahavior of emax is shown in Fig. 5 as a function 
of Pvac/p- Galactic clusters with Pvac/p ^ 0.1 may have a vanishing mean velocity for e^ax — ?■ 1 in the 
oblate case and e^ax ^ 0.92 in the prolate case. If the density is smaller, a vanishing mean velocity can 
be reached for non so fiat objects. Spherical configurations have (v^) for p 2/9vac- 
Since the mean velocity is proportional to the temperature, p = M/V , Eq. represents also a 

mass-temperature relation. Hence the results for the mean velocity squared are also applicable to the 
temperature of the configuration. For instance, we can write Eq. ^3 for a cosmological structure, 
say a galactic cluster, by writing its density as resulting from a perturbation Sp from the background 
density of the universe pYj(t) as pgc = Pb(^)(l +<^(^))) where S(t) = Sp(t)/pYj(t). Equation allows us 
to determine the temperature of the cluster at a given cosmic time t as 



T 



Tie) 



lOfc, 



\-A[e) 



l-fiv 



(71) 



e^ arcsin e 



Prolate, 
Oblate, 



where a[t) is the scale factor and p is the mass of average components of the cluster. Equation ^3 
assumes a fiat universe Slmatter + f^vac = 1- This result is a generalization of a result derived in ^3 
valid for spherical geometry. In that case one recovers the typical mass temperature relation T cx M^/^ 
mantaining p constant. Although this has the same dependence as in equation the meaning is 

different since ^3 for A = ^ computes the temperature of a certain galactic cluster at some redshift 
given its mass while equation is associated to a reversible process where a configuration passes from 
T = to some final through states of virial equilibrium mantaining a constant mass. 



5 Small oscillations in the Newton-Hooke spacetime 

The stability condition of Newtonian configurations against oscillations can be also derived from the 
second order virial equation by expanding the periodic Lagrangian perturbations £(r,t) = £(r)e''^* with 
oscillation frequency cj. For simplicity, we will consider a rotational configuration without internal 
motions. By assuming adiabatic perturbations, the variational form of the second order virial equation 
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El is written as (see ^ElEI 

1 A o 16 
~ 2 



+ l^^kj^i^ - 3r)p [{n ■i)in-r)- fi2(r . £)] d 



where F = (9 In _P/(9 Inp) ^ is the adiabatic index governing the adiabatic perturbations. For polytropic 
e.o.s we have 7 = F, of course. The Fagrangian displacement is constrained with the boundary conditions 
i), £ = at r = and ii) £ must be finite at the surface. The boundary conditions are required to solve 
the equation of motion for £ resulting from perturbing Euler's equation (the Sturm-Fiouville eigenvalue 
equation). 

Fet us consider radial oscillations which are adequate for spherical systems and for small deviation from 
spherical symmetry. By assuming a general trial function for the Fagrangian displacement £ = /(r)r 
satisfying the boundary conditions and taking the trace in Eq. we can solve for frequency of the 

oscillation about equilibrium 



2 



3^3; 3 \ Q 



(73) 



where we have assumed F = constant throughout the configuration and we have defined in analogy to 
Eq. ^3 the quantities 

\W\ = -Q^\W\, I = Qi, A = —^. (74) 
' ' 2^ ' ^ ' 3 \W\ 

In Eq. ^3 g is a parameter with units of density, and 

f d$^ If f 

W=- / p{r)f{r)-—d\, R=-nlJ p{r)f{r)r-'[P -SijXiXj]d\, 1= / p{r)f{r)rd\, 
Jv cl*" ^ Jv Jv 

(75) 

together with B = The critical adiabatic index is written as 



Tent = ^ (1 - 25)-' 



l + iA^--B 
2 g 2 



(76) 



such that for F < Fcrit instability sets in and the system becomes unstable while perturbed. From 
Eq. 131, we see that through the inclusion of the cosmological constant we (the system) are forced to 
choose a bigger adiabatic index. In the simplest case, when /(r) = r, we have W — ?> W, / — ?> I, 5 — ?> /3. 
Furthermore, for p =constant, we get A ^ A and g = p. The stability condition then becomes 



r> ^(1-2/3)-' 



2 p 2 



(77) 



The spherical symmetry which we assumed compels us to write the ratio /3 for low eccentricities or to 
parametrize it in terms of the total angular momentum of the configuration. In the first approximation, 
at low eccentricities the ratio /3 given in takes the form 

/3^^Q-^)e^+0(e^). (78) 

This equation is useful if we want to calculate a small, but peculiar effect. Insisting that F is very close 
to 4/3 we can convert the stability condition F > Fcrit into a condition on eccentricity, namely 

e > e^,„ = 3.35 (^1 + 5^ - 3^^ + ' (^9) 

where e = | — F measures a small departure from the critical value 4/3. For the specific case e = 0, we 
conclude that the eccentricity must be such that emin < e ^ 1 with 

e^,„ = 3.8 (^^y^' , (80) 
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Figure 6: Critical adiabatic index as a function of /? for a homogeneous spherical configuration for different values of 

C = 2pvac/p 

which is clearly valid only for large densities, say for p > lO^pvac- This imposes a range for the 
eccentricity ecrit < e ^ 1 for stability under radial perturbations. 

For completely spherical configurations with constant density, we can write fi = Vt^^^/ATrp. Figure 
shows the behavior of the critical adiabatic index as a function of fi for different values of C = 'Ip^^d P- 
The largest deviation is for low densities as expected. 
For non rotating configurations, the stability criteria is simply given by 



1 Pvac 

2 e 



A 



Stt p{r)P f{r)dT 



(81) 



If we work with a constant density profile, the choice of another trial function satisfying the boundary 
conditions mentioned before will not play a role in the above expressions since A ^ A = 2 for any f{r). 

Related results regarding the adiabatic index in cosmologies with non-zero A have been obtained in the 
relativistic frame work in |^|. 

6 Conclusions 



In the present work we investigated in detail the astrophysical relevance of the cosmological constant for 
equilibrium configurations. Using the tensor and scalar virial equations and the Lane-Emden equation 
we could show that many astrophysical facets get modified by A. The second aspect concerns the fact 
that A introduces new relevant scales (these scales would be zero if A = 0) like the maximal virial 
volume defined by the maximal virial radius and the maximal extension of bound orbits given in 
Eq. E3- 

It is often assumed that superclusters with densities ~ pcrit are not in equilibrium. With the inequality 
^3 we have a precise tool to quantify this statement. Indeed, the pancake structure of the superclusters 
lead us to the conclusionthat they are even far away from the equilibrium state due to the factor A 
which grows with the object's fiatness. On the other hand, relatively low density objects can still reach 
equilibrium (even if ^ = 2 for the spherical case, the density of the object has to be only twice as large 
as Pvac, i.e, 1.4/9crit). There is nothing which could, in pronciple, prevent a relatively low density object 
to be in equilibrium This is not only a result obtained from the virial equations with A, but follows also 
from hydrostatic equilibrium with the inclusion of A It is then natural to put forward the question 
how the shape of such objects affect their properties as compared to the case without A. For instance, 
if p/ Pvac = 10 and e ~ 0.9, the angular velocity is reduced approximately 30% with respect to its value 
with A = (Slrot = O.TSlo). Other effects discussed in the text are the mean velocity (see Fig. 5) and 
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the adiabatic index (see Fig. 6). To be specific let us takje an example of a very low density object, 
p = 2.566/9crit- This object is still virialized as long as its eccentricity is e < 0.5. Approaching e 0.5 the 
mean velocity (temperature) goes to zero until finally beyond 0.5 the objects ceases to be in equilibrium. 

To highlight some quantities which get affected let us mention the polytropic index n and the 
angular velocity. Low-density objects with p ~ 20/9crit an index n > 1.5 becomes unphysical as the the 
conglomeration of matter does not have a definite finite radius. Taken together with the results of section 
five on stability against small oscillations this implies that the adiabatic index 7 from the equation of 
state P cc p"' gets restricted for these densities in a rather strong way, namely by 1.33 < 7 < 1.66. 
For higher densities the effect grows with n. For instance, for p ~ lO'^pcrit, « = 5 is an unacceptable 
solution due to the lack of a finite extension. For non-spherical configurations we could show an effect 
of A on angular velocities, the mean velocity of the component of the astrophysical object and on the 
critical mass, again for densities one and two orders of magnitude above the critical one. We think that 
the work can be generalized in various ways. For instance, to generalize the virial equations to be able 
to differentiate between different models of dark energy mentioned at the beginning. Another aspect is 
to look into the Lane-Emden equation for non-spherical geometries ^3^3 • Here we concentrated on 
the equilibrium condition of already virialized matter. The collapse of matter in an accelerated universe 
with dark energy has been considered in 



Appendix 

Here we briefiy show the solution for the cubic equation . Let us write that expression as 

t/^+ra + g = 0, p=10?yri>0, g = -3r,ri<0, (82) 

corresponding to a positive cosmological constant. Associated with this expression one defines a dis- 
criminant and a quantity R defined as 



^=^/ + |g'>0, R = sgn{q)^^\p\<0. (83) 
Since p > 0, the roots of are given in term of the auxiliary angle (f) defined as 

smh,/-=^, (84) 
so that the only real solution for is written as 

(85) 



y = -2i?sinh ( ^ j = ^^\p\smh 



larcsmh (^) 
3 \2R3J 



The other two roots are complex numbers whose real part are negative. If p < and q > the roots 
depend whether D is bigger or smaller than zero. If _D < 0, the three roots are real and negative, while 
D > yields one real negative root and two complex roots which real parts are given as 



— \p\ cosh 
3 



larccosh (^) 
3 \2R3J 



(86) 
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